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CN ■ Primordial gas in protogalactic dark matter (DM) halos with virial temper- 

atures Tvii- > lO^K begins to cool and condense via atomic hydrogen. Provided 
this gas is irradiated by a strong ultraviolet (UV) flux and remains free of H2 
^ ■ and other molecules, it has been proposed that the halo with Tvir ~ lO'^K may 

^ . avoid fragmentation, and lead to the rapid formation of a supermassive black 

P^i hole (SMBH) as massive as M ^ 10^ - 10^ M©. This "head-start" would help 

Q . explain the presence of SMBHs with inferred masses of several x IO^Mq, powering 

-)— > . the bright quasars discovered in the Sloan Digital Sky Survey at redshift -z > 6. 

c^ . However, high-redshift DM halos with T^ir ~ lO^K are likely already enriched 

with at least trace amounts of metals and dust produced by prior star-formation 
^ , in their progenitors. Here we study the thermal and chemical evolution of low- 

metallicity gas exposed to extremely strong UV radiation flelds. Our results, 
obtained in one-zone models, suggest that gas fragmentation is inevitable above 
^ , a critical metallicity, whose value is between Z^r ~ 3 x 10~^Zq (in the absence 

^ ! of dust) and as low as Zcv ~ 5 x 10~^Zq (with a dust-to-gas mass ratio of about 

00 ! 0.01Z/Zq). We propose that when the metallicity exceeds these critical values, 

dense clusters of low-mass stars may form at the halo nucleus. Relatively mas- 



^ 



sive stars in such a cluster can then rapidly coalesce into a single more massive 
object, which may produce an intermediate-mass BH remnant with a mass up 



.<^ ; to M < 10^ - 10^ Mq. 
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Introduction 



The discovery of bright quasars at redshifts z > 6 in the Sloan Digital Sky Survey 
(SDSS) implies that BHs as massive as several xlO^M© were already assembled when the 
age of the universe was less than ^ 1 Gyr (see the recent review by Fan 2006). The BH 
masses are inferred from the quasars' luminosities, assuming these sources shine near their 
Eddington limit. Strong gravitational lensing or beaming could, in principle, mean that the 
inferred BH masses are overestimated; however, there is no obvious sign of either effect in 
the images and spectra of these quasars (Willott et al. 2003; Richards et al. 2004). Indeed, 
their relatively "normal" line-to-continuum ratio, consistent with those in lower-redshift 
quasars, makes it unlikely that the apparent flux of these sources was significantly boosted 
by beaming (Haiman & Cen 2002). Likewise, the lack of a second detectable image on Hubble 
Space Telescope images (Richards et al. 2004) essentially rules out the hypothesis that most 
of the sources experienced strong magnification by lensing (Comerford et al. 2002; Keeton 
et al. 2005). 

Relatively little time is available for the growth of severalxlO^ Mq SMBHs prior to 
z ^ 6, and their seed BHs must be present as early as 2; ~ 10 (e.g. Haiman & Loeb 
2001). As the SMBHs grow from high-redshift seed BHs by accretion, they are expected 
to encounter frequent mergers. A coalescing BH binary experiences a strong recoil due to 
gravitational waves (GWs) emitted during the final stages of their merger. The typical recoil 
speed is expected to be f recoil ^ 100 km s~^ (and may be as large as 4, 000 km s~^ for special 
BH spin configurations; see, e.g. Campanelli et al. 2007 and references therein), significantly 
exceeding the escape velocity (< 10 km s~^) from typical DM halos that exist at z ~ 10. As 
a result, SMBHs are often ejected from their host halos at high redshift. The repeated loss of 
the growing seeds makes it especially challenging to account for the severalxlO^ Mq SMBHs 
at 2; > 6 without at least a brief phase of super-Eddington accretion, or some equivalent 
"head-start" (Haiman 2004; Yoo & Miralda-Escude 2005; Shapiro 2005; Volonteri & Rees 
2006). 

There have been several recent proposals that such a "head-start" may occur in metal- 
free gas in high-redshift DM halos with virial temperatures exceeding T^u > lO^K, leading 
to the rapid formation of SMBHs with a mass of M ^ 10^ — 10^ Mq. As primordial gas falls 
into these halos, it initially cools via the emission of hydrogen Lya photons. Provided the gas 
is free of H2 molecules, its temperature will remain near T^^,. ~ lO^K. Bromm & Loeb (2003, 
hereafter BL03) performed hydrodynamical simulations of a metal- and H2-free halo, with 
a mass of ~ IO^Mq collapsing at 2; ~ 10, corresponding to a 2a Gaussian overdensity and to 
Tvir ~ 10^ K. Under these conditions, which may apply to some dwarf galaxies collapsing close 
to the epoch of reionization, the primordial gas is marginally able to collapse and remains 
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nearly isothermal. BL03 found that during the evolution, fragmentation of the gas cloud 
is very inefficient, leading at most to binary formation even with some degree of rotation. 
Thus, a super-massive star is expected to form, and evolve into a SMBH with a mass as 
high as M ^ 10^ - 10^ Mq. Oh & Haiman (2002) and Lodato & Natarajan (2006) have also 
showed that if H2 formation is inhibited, a primordial-gas disk is stable to fragmentation and 
a single massive object is formed in accordance with BLOS's conclusion. Volonteri & Rees 
(2005) arrived at similar conclusions, by considering Bondi accretion onto a stellar seed BH, 
which can significantly exceed the Eddington rate at the gas density and temperature in a 
similar halo. Finally, Begelman et al. (2006) and Spaans & Silk (2006) proposed different 
mechanisms to form similarly massive BHs by the direct collapse of primordial, atomic gas. 
For reference, we note that the total (DM+gas) mass of halos with Tvir = (1 — 5) x lO'^K at 
2 = 10 is Mtot ^ 10^-9 M0, so that such SMBHs would represent ^ 0.2 - 20% of the gas 
mass in these halos. We also note that in the WMAP5 cosmology, the age of the universe 
at 2; = 10 and z = 6.5 is ~ 0.5Gyr and ~ 0.9 Gyr, respectively. At the e-folding time-scale 
of 4 X 10^ years (assuming Eddington accretion, and a radiative efficiency of 10%; see, e.g., 
Haiman & Loeb 2001), a seed BH of M ^ IO^Mq at 2; ~ 10 could easily grow to a super 
massive BH of M ^ 2 x IO^Mq at 2; ~ 6.5, if fed uninterruptedly. 

A crucial assumption in all of the above proposals is that H2 molecules cannot form as 
the gas cools and condenses in the DM halo. This assumption can be justified in the presence 
of a sufficiently strong far ultraviolet (FUV) radiation, so that molecular hydrogen (or the 
intermediary H~ necessary to form H2) is photodissociated. The relevant criterion is that 
the photodissociation timescale is shorter than the H2-formation timescale; since generically, 
tdiss oc J and tfomi oc p, the condition tdiss = ^form yields a critical ffux J oc p. In DM halos 
with Tvir ^ lO'^K, whose gas can not cool in the absence of H2, the densities remain low 
and H2 can be dissociated even when background ffux is as low as J21 ~ 10~^ (e.g. Haiman, 
Rees & Loeb 1997; Mesinger et al. 2006; here J21 is the flux just below 13.6eV, in the usual 
units of 10~^^erg cm~^ sr^^ s~^ Hz~^). However, if a gas cloud is massive enough and has a 
virial temperature higher than ^ 8000 K, it is able to cool and start its collapse via atomic 
hydrogen Ly-a cooling. Even if the FUV ffeld is initially above the critical value, molecular 
hydrogen can form, and dominate the gas cooling at a later stage during the collapse (Oh & 
Haiman 2002); the H2-formation rate is furthermore strongly boosted by the large out-of- 
equilibrium abundance of free electrons in the collisionally ionized gas in these halos (Shapiro 
& Kang 1987; Susa et al. 1998; Oh & Haiman 2002). The critical ffux required to keep the 
gas H2-free as it collapses by several orders of magnitude therefore increases signiffcantly; 
for halos with T^ir ~ lO'^K the value has been found to be J21 ~ 10^ — 10^, depending on the 
assumed spectral shape (Omukai 2001, hereafter O2001; BL03). In halos exposed to such 
extremely intense UV ffelds, the gas cloud is still able to collapse only via atomic hydrogen 
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line cooling, namely Ly a and H free-bound (f-b) emission (O2001). 

One possible source of such an intense UV field is the intergalactic UV background 
just before the epoch of cosmic reionization (BL03). The ionizing photon flux J21 can be 
evaluated from the number density of hydrogen atoms in the intergalactic medium (IGM) 
and the average number of photons needed to ionize a hydrogen atom A^^, which, in general, 
is > 1, owing to recombinations in an inhomogeneous IGM. Using the escape fraction of 
ionizing radiation /esc, the flux J^ just below the Lyman limit is given by 

'' /esc /esc47r run \10J Vo.Oi; V H 

where Yh = 0.76 is the mass fraction of hydrogen, mn is the proton mass, and pb is the baryon 
density (assumed here to correspond to fife/i^ = 0.023; Dunkley et al. 2008). Equation [1] 
shows that J^ can approach the critical value at z > 10, provided /esc is small, and N^ is 
large; fesc/Nj < 10^'^. Although the value of /esc is quite uncertain, in low-mass minihalos, 
the expectation is /esc ~ 1, as these halos are easily self-ionized, and most of their ionizing 
radiation escapes into the intergalactic medium (Kitayama et al. 2004; Whalen, Abel & 
Norman 2004). Observations of nearby star forming galaxies indicate lower values /esc < 0.1 
(Leitherer et al. 1995, Inoue et al. 2005); nevertheless, it appears unlikely that the mean 
cosmic background will reach the critical values. A few halos that have close and bright 
neighbors may still see a sufficiently increased flux (Dijkstra et al. 2008). Alternatively, the 
critical value could be established by sources internal to the halo, e.g. by a vigorous phase 
of starburst (Omukai & Yoshii 2003) or by the accreting seed BH itself (Volonteri & Rees 
2005). 

So far, the process of supermassive star (and SMBH) formation in dwarf galaxies in 
the presence of a strong UV background has been investigated only under the hypothesis 
that the gas is metal-free. However, the T^ij- ~ lO^K, M ~ 10^ Mq halos forming close to 
the epoch of reionization are built from lower-mass progenitors that had collapsed earlier. 
Many, and perhaps all of these halos should therefore be enriched with at least some trace 
amount of metals. Furthermore, is it unlikely that the strong critical FUV flux could be 
generated and maintained without efficient star formation at higher redshifts (if J21 ^ 10^ 
was produced by accreting BHs, this would signiflcantly overpredict the present-day soft 
X-ray background; Dijkstra et al. 2004). It is well-known that adding metals and dust into 
a primordial gas, even at trace amounts as low as Z ~ 1O~^Z0, can signiflcantly affect its 
cooling properties (Schneider et al. 2003; Omukai et al. 2005; Schneider et al. 2006). Two 
independent hydrodynamic simulations (Tsuribe & Omukai 2006; Clark, Glover, & Klessen 
2008) recently studied the fragmentation of metal-enriched collapsing protogalactic clouds, 
in the absence of an external FUV fleld, and found efficient fragmentation for Z > IO^^Zq. 
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In the present paper, our goal is to answer the following question: can cooling and 
fragmentation be avoided in metal-enriched T^it- > l^i^K halos, irradiated by a strong FUV 
flux? If so, this would suggest that supermassive black holes may form, similar to the metal- 
free case, in the more likely case of metal-enriched high-redshift protogalaxies. To investigate 
this possibility, we here study the thermal and chemical evolution of low-metallicity gas, 
exposed to extremely strong UV radiation fields. We will evaluate the critical metallicity, 
above which fragmentation becomes unavoidable in the presence of a strong FUV flux. 

In § El we describe our one-zone modeling procedure. Our results are presented and 
discussed in §[3], first for the metal-free (§ 13. ip . and then for the metal-enriched case (§ 13. 2p . 
The fragmentation and subsequent evolution of the metal-enriched clouds are then discussed 
in § 13.31 and 13.41 respectively. In § IH we summarize our results and offer our conclusions. 



2. Model 

2.1. Basics 

We use the one-zone model described in Omukai (2001) to follow the gravitational 
collapse of gas clouds. The model includes a detailed description of gas-phase chemistry and 
radiative processes, and the effect of dark matter on the dynamics in a simplified fashion. In 
addition, in the present version of the model we have implemented the contribution of metal 
lines and dust to gas cooling. 

In what follows, all physical quantities are evaluated at the center of the cloud. The gas 
density increases as 

^Pgas _ Pgas /2\ 

dt tcoi 
where the collapse timescale, tcoi, is taken to be equal to the free-fall time. 



and p is the sum of the gas and dark matter density. The dark matter density follows the 
evolution of a top-hat overdensity, 

PDM- — (^Y^^^J ^DMPcrit (4) 

with 

1+2 = (1 + Zta)f j (5) 
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(e.g., Chapter 8.2 of Padmanabhan 1993), where the turn-around and the viriahzation cor- 
respond to ^ = TT and 27T, respectively. Although, strictly speaking, this is correct only in 
the Einstein-de Sitter universe [Qq = 1), it does not cause a significant error in the high-z 
universe {z > 10) we consider. 

The initial epoch of calculation is taken at the turn-around at redshift Zta = 17. From 
equation [5l the viriahzation and turn-around redshifts have the relation 1 + z^-^j. = 2~^/^(l + 
Zta); thus 2;vir — 10. In our calculation, the dark matter density is kept constant after reaching 
its viriahzation value 8pDM(-2^ta)- The initial values of the gas number density, temperature, 
ionization degree, and H2 fraction have been assumed to be wh = 4.5 x 10~'^cm~^, T = 21 K, 
y{e) = 3.7 X 10~^ and y{B.2) = 2 x 10~®, respectively, to refiect conditions at the turn-around 
at Zta = 17. Some runs with initial temperature ten times higher (210K) are also performed 
to confirm independence of our main results from the initial temperature. The cosmological 
parameters are Qbm = 0.24, Q\, = 0.04, and h = 0.7. 

Our calculation does not include the viriahzation shock. Owing to fast cooling by Lya 
emission, the central region whose evolution we intend to follow does not experience the 
viriahzation shock in the spherically symmetric case (Birnboim & Dekel 2003). In more real- 
istic calculations, the outer regions can experience shocks and the temperature and electron 
fraction become higher than in our case. In addition, recent numerical calculations (e.g. 
Keres et al. 2005) show that low-mass galaxies, especially at high-redshifts, obtain their gas 
through accretion predominantly along the large-scale filaments. Three-dimensional effects 
such as asymmetric accretion might affect the evolution at low densities. However, since we 
are considering halos with Tvir — lO'^K, which can marginally collapse by Lya cooling, the 
shock is not strong: the temperature increase is modest and the electron fraction reaches 
at most < 10^^ (see Figures 5a and 5c in BL03). This additional electrons alter the early 
evolution for the J = case. However, in the irradiated clouds, where H2 formation is 
suppressed, during the collapse by the Lya cooling recombination proceeds until the free 
electron fraction reaches Xe — 1.2 x 10~^n^ , the value set by the balance between the 
recombination and the collapse time tree ~ tcoi at 8000K. Thus, our results for molecule 
formation and cooling are hardly affected. 

We adopt tg as the collapse time scale just because it has been widely used in other 
studies (e.g., Palla et al. 1983). Note that the free-fall time ([3]) is the time for density of an 
initially static cloud to reach infinity, while the dynamical timescale tcoi = p/{dp/dt) in the 
free-fall collapse is 

in the limit where the density has become sufficiently larger than the initial value. Thus the 
rate we adopted ([3]) is 37r/2 = 4.7 times slower collapse than the genuine free-fall one. In 
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fact, pressure gradients oppose gravity and the collapse becomes slower than the free-fall one 
within a factor of a few (e.g., Foster & Chevalier 1993). Adoption of t^ as the e-folding time 
for density increase mimics the pressure effect. The assumption of nearly free-fall collapse 
is invalidated, and the collapse is slowed down, once the cloud becomes optically thick to 
continuum radiation. However, our result on the thermal evolution is not altered: with little 
radiative cooling, the temperature is now determined by the adiabatic compression and the 
chemical cooling by dissociation and ionization, both of which are independent of the collapse 
timescale. Moreover, the evolution after the cloud becomes optically thick is not relevant to 
our argument on fragmentation, which occurs at much lower density, in the optically thin 
regime. 

The overall size of the collapsing gas cloud (or of the roughly uniform density central 
region) determines its optical depth, and is therefore important for its thermal evolution. 
Here we assume the size equals the Jeans length. 



where Tgas is the gas temperature, fi is the mean molecular weight. Similarly, its mass is 
given by the Jeans mass 

Mj = PgasAJ^. (8) 

Specifically, we assume that the radius of the cloud is Re = Aj/2 and the optical depth is 

T^ = KyRc = '^i^ ( Y ) • (9) 

In addition to dust absorption (see § l2.2p . we include the following primordial-gas processes 
as sources of the opacity Kj^ (Table 1 of O2001): the bound-free absorption of H, He, H~, 
Hj', free-free absorption of H~, H, collision-induced absorption of H2-H2 and H2-He, Rayleigh 
scattering of H, and Thomson scattering of electrons. 

The temperature evolution is followed by solving the energy equation: 

de d f 1 \ Anet 



^ = _p^ _ _— , (10) 

dt dt VPgas/ Pgas 

where e is the internal energy per unit mass 

e = ^— ^ (11) 

7ad - 1 yumn 

p is the pressure, 7ad is the adiabatic exponent, and Anet is the net cooling rate per unit 
volume. In addition to cooling and heating processes for the primordial gas, which include 



continuum emission, as well as emission by H and H2 lines, and chemical heating/cooling, 
the net cooling rate includes emission by C and O fine-structure lines Amctah by dust grains 
Agr, and heating by photoelectric emission of dust grains Fpe. Cooling by fine-structure fines 
of [CII] and [01] is included as in Omukai (2000). Dust processes are described below in 

Primordial-gas chemical reactions are solved for the nine species of H, H2, e, H"*", H^, 
H~, He, He"*", and He"'""'". We do not explicitly include the chemical reactions involving metals. 
Instead, all the carbon and oxygen is assumed to be in the form of CH and 01, respectively. 
Having a lower ionization energy (11.26 eV) than hydrogen, carbon remains in the form 
of CH in the atomic medium owing to photoionization by the background radiation. We 
maintained this assumption even in J = runs, although carbon is expected to recombine 
and become neutral in these cases. The cooling rates by CII and CI fine-structure lines 
are within a factor of ~ 2 difference for T > 30K, and therefore this assumption does 
not significantly affect the results. On the other hand, the ionization potential of oxygen 
(13.61 eV) is very similar to that of hydrogen (13.60 eV) and the charge exchange reaction 

0+ + H^H+ + 0, (12) 

keeps its ionization degree equal to that of hydrogen. In fact, the coefficient of the rightward 
reaction being 6.8 x 10~^°cm^/sec, these reactions reach equilibrium only in ~ SOra^^yr. 

In a cold (< a few lOOK) and dense (> lO^^^cm^^) environment, molecular coolants 
such as CO and H2O may become important (Omukai et al. 2005). Since we are interested 
here in metal effects on warm (> a few lOOOK) atomic clouds, we neglect the contribution to 
cooling of metals in molecules. This simplification does not affect the early evolution of gas 
clouds, when the effects of metals induce a deviation from the primordial evolutionary track 
at several lOOOK. It is true that it may alter the predicted thermal behavior at later stages, 
when the gas has cooled significantly (< lOOOK). However, even in such cold environments, 
the error in the temperature caused by neglect of metal molecular coolants is very small (see 
Figure 10 of Omukai et al. 2005) and the thermal evolution is well reproduced when only 
dust processes and fine-structure line cooling of C and O are considered. 



2.2. Dust Processes 

Dust in the local interstellar medium (ISM) originates mainly from the asymptotic 
giant-branch (AGB) stars, whose age is > 1 Gyr, longer than the Hubble time at z > 6. At 
higher redshifts, supernovae (SNe) are considered to be the major dust factories. Indeed, 
the observed extinction law of high-2; quasars and gamma-ray bursts can be well reproduced 
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by this scenario (Maiolino et al. 2004, Stratta et al. 2007). Dust grains produced in SN 
ejecta are more effective in cooling and H2 formation because of tlieir smaller size and larger 
area per unit mass (Schneider et al. 2006). However, their composition and size distribution 
are still affected by many uncertainties, such as the degree of mixing in the ejecta and the 
efficiency of grain condensation and their destruction by the reverse shock (Nozawa et al. 
2007, Bianchi & Schneider 2007). 

To be conservative, in this work the properties of dust, such as grain composition and 
size distribution, are assumed to be similar to those in the solar neighborhood and its amount 
is reduced in proportion to the assumed metallicity of the gas clouds. Specifically, we adopt 
the dust opacity model developed by Semenov et al. (2003). This model partly follows the 
scheme proposed by Pollack et al. (1994), which was used in Omukai et al. (2005), assuming 
the same dust composition, size distribution and evaporation temperatures, but uses a new 
set of dust optical constants. Overall, the opacity curves of the two models are in good 
agreement, the largest difference being at most a factor of two (see Semenov et al. 2003 
for a thorough discussion). The main dust constituents include amorphous pyroxene ([Fe, 
MgjSiOs), olivine ([Fe, Mg]2Si04), volatile and refractory organics, amorphous water ice, 
troilite (FeS) and iron. The grains are assumed to follow a size distribution modified from 
that by Mathis, Rumpl, & Nordsieck (1977) with the inclusion of large (0.5 - 5)/im grains. 

At each density and gas temperature, the dust is assumed to be in thermal equilibrium, 
and its temperature Tgr, which is followed separately from the gas temperature, is determined 
by the energy balance equation 



47r / Ks,,uB,y{Tgj.)dU = Agas^duSt + 4:71 Ks.^yX^dv. 



(13) 



Here Agas^dust is the energy transfer rate per unit mass from gas to dust due to gas-dust 
collisions, which we take from Hollenbach & McKee (1979), Ka.jy is the absorption opacity 
of dust, and J™ is the mean intensity of the radiation field inside the cloud. Note that 
Agas^dust also represents the net cooling rate of the gas, caused by the presence of dust 
grains at temperature Tgr. We model the external radiation field assuming a diluted thermal 
spectrum (i.e. a blackbody spectrum, scaled by an overall constant representing a mean 
geometrical dilution). Its shape is then fully described by only two free parameters, J21, the 
mean intensity at the Lyman limit (z/h) and T^,, the color temperature, 

JT = J2il0-'M5,(T,)/5,JT,)] erg cm-^sr-^s-^Hz-^ (14) 

In the following, we will consider two possible values for the radiation color temperature, T^, = 
10^ K and 10^ K, representing "standard" Population II stars and very massive Population 
HI stars, respectively. Given the mean intensity of the external radiation field J^^, the field 
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inside the gas cloud is obtained as (see O2001), 

7CX _|_ t ™ c 
Tin _ '^v ^ ^u-^i"~'a.,u ,-. p.^ 

■^^ " l + ^.x. ' ^^^' 

where 1 — ^,y is the single-scattering albedo, S'a,,^ is the source function, 

Xj, = max [r^,rj, (16) 

and the optical depth r^ includes both dust and gas opacities. 

The gas is heated by photoelectrons ejected from dust grains after absorption of FUV 
photons. Following Bakes and Tielens (1994), we compute the photoelectric heating rate as, 

r(f) = Tpe-Ape (17) 

= [lO-^^eGo^H - 4.65 x 10-^''T'^-''\GoT^/yn{e))''n{e)nn]Z/ZQ (18) 

where 

4.9 X 10-^ 3.7 X 10-^(T/10-^)°-^ 

^ ~ [1 + 4 X 10-3(G'oTi/7n(e))0-73] ^ [1 + 2 x IQ-^GoT^ynie))] ' ^ ' 

and (3 = 0.735/T°'°^^. The Habing parameter Go is defined as 

Go = 47r / J^du /l.6 x 10-^ (20) 

where the integral is over the FUV radiation in the range (5.12-13.6) eV. In the above formula 
(eq. [18]), small (< 15A) grains contribute about half of the heating. 

Molecular hydrogen efficiently forms on dust grains by depositing the chemical binding 
energy on the grain surface. Following Tielens and Hollenbach (1985), we compute the H2 
formation coefficient as, 

, ^ 6.0 X IQ-'-^jT/SOOKy/^UZ/ZQ 

^' 1 + 4.0X 10-2(T + Tg,)i/2 + 2.0x 10-3T + 8.0X 10-6T2 ^ ^ 

where 

/a = [1 + exp(7.5 X 10^(1/75 - l/T^r))]'^. (22) 



3. Results 

In what follows, we will first discuss the results obtained for the thermal evolution of 
metal-free gas clouds, and then describe the effects induced by the presence of metals and 
dust grains. 
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3.1. Metal-free Clouds 



The thermal evolution of metal-free clouds irradiated by a FUV radiation background 
is expected to change with radiation temperature T^ and intensity J2i- The models with 
a radiation temperature of T* = 10^ K (10^ K) are shown in Figured] ([2l respectively) for 
different values of intensity J21. 

Initially, i.e. at low densities, the temperature increases adiabatically, because there is 
not enough H2 to activate cooling. In the no radiation case, when the density is ~ lcm~^ 
and the temperature is ~ 1000 K, sufficient H2 is formed and, as a result, the temperature 
decreases. It is to be noted that the relatively low temperature where this condition is met 
does not contradict previous results (BL03). In fact, the predicted temperature of each fluid 
element in the simulation of BL03 shows a large scatter at low densities. This scatter reflects 
the radial temperature gradient, and the central value, which we calculate here, corresponds 
to the lower boundary of the scattered points and it is in agreement with our result. We 
expect that the central temperature of the gas cloud does not reach the virial temperature 
of the host dark matter halo since the innermost region starts to cool and collapse during 
the adiabatic compression and does not experience the virialization shock. 

As the external radiation intensity J21 increases, the onset of H2 cooling is delayed 
because higher densities and temperatures are required for H2 formation to compensate for 
the photodissociation. If the UV intensity is below a threshold value, J2i,thr, which we flnd 
to be in the range 10^ - 10^ for T* = 10^ K and (1 - 3) x 10^ for T, = 10^ K, there is always 
a density at which II2 cooling starts to become effective. The temperature then decreases 
and eventually reaches the no-radiation evolutionary track, along which it evolves thereafter. 
On the other hand, if the radiation is stronger than the threshold value, H2 cooling never 
becomes important. In this case, atomic hydrogen cooling by H excitation (for < lO^cm"^) 
and H~ free-bound (f-b) emission (for > lO^cm"^), are the main cooling channels (see Figure 

ED. 

In Figure dl runs with higher initial temperature (210 K) are also shown (dotted lines). 
During the initial adiabatic phase, the temperature at a given density is proportional to its 
initial value, and thus higher in runs with higher initial temperature. However, after the 
onset of efficient radiative cooling, these initially different thermal evolutionary tracks soon 
converge. At higher densities, the results are independent of the initial temperature (see 
Figure [T]) . 

As it can be inferred from Figs. [1] and [21 we flnd that the threshold value, J2i,thr, is 
lower for a radiation temperature of T* = 10^ K than for T* = 10^ K. Thus, for comparable 
radiation intensities, J21, the lower T^ radiation has a stronger impact on the cloud evolution. 



-12- 



To understand why this is the case, in Figure H] we show the H2 and H^ photodissociation 
rates, for the same intensity J21 = 1. The dilution factor W, defined by Ji^ = WB^(T^), 
which was used in Omukai & Yoshii (2003), is also shown for reference. As the figure shows, 
the H~ photodissociation rate decreases steeply with T^,, while the H2 photodissociation rate 
remains nearly constant. The H2 and H~ photodissociation rate coefficients are 

fcH^ph = 1.4 X 10V^(12.4eV) (23) 

in the unattenuated case and 

^H-ph = / —-^a^du (24) 

where the lower limit on the latter integral, 0.755 eV, is the threshold energy for H~ pho- 
todissociation. Since the radiation field is normalized at the Lyman limit (13.6 eV), ^Haph is 
not sensitive to T*, whereas /cn-ph, which refiects the radiation field above 0.755 eV, depends 
significantly on the adopted T^. H2 formation proceeds via a two step process (H~ channel), 

H + e^H-+7, (25) 

and 

H- + H ^ H2 + e. (26) 

If H~ is photodissociated, it can not activate the second step fl2B]) . This is the reason why a 
lower T* radiation field leads to a less efficient H2 formation rate and to a lower H2 fractional 
abundance than a higher T^, radiation field. 

Cooling via H^ f-b emission occurs through the radiative association reaction (!25|) . 
Subsequently, H~ is collisionally dissociated through the reaction 

H- + H ^ 2H + e, (27) 

and the whole process results in a net cooling by the emitted photon. Thus, as long as the 
collisional dissociation rate exceeds the photodissociation rate, the H~ f-b emission is hardly 
quenched, even in the presence of a strong FUV radiation field. Due to the small opacity, 
H~ f-b cooling becomes important only at high densities (> lO^cm"^) and temperatures (~ 
several 10^ K), where the above condition is always satisfied. Then, H~ f-b cooling is not 
affected by photodissociation. 

As shown in Figure [31 radiative cooling and compressional heating rates almost balances 
for a wide range of densities. Namely, in equation (llOp . the two terms in the right hand side 
are almost cancels and the left hand side is much smaller than those terms. Suppose that 
the radiative cooling rate per unit volume depends on the density and temperature as 

Arad/Pgas OC p^^^'T^ . (28) 



-13- 

Since the compressional heating rate 
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oc p'£t, (29) 



'2'' Pgas Pgas''col 

the thermal balance of those two terms results in the following temperature evolution: 

Tocpgas^. (30) 

Both the H line and H~ f-b emissions are very sensitive to temperature, and thus /3 > 1. For 
those collisional processes, a = 2 for fixed chemical abundances. With chemical evolution, 
it deviates from 2, but remains > 3/2. Therefore, the exponent in equation (13U]1 is negative 
for the atomic-cooling track as long as the cloud is optically thin: the temperature decreases 
with density as observed in Figures [1] and [2l On the other hand, on the molecular-cooling 
track, a = 1 for densities higher than the critical value for the LTE. Thus, the temperature 
increases with density for rin ^ lO^cm"^. 

The existence of a threshold UV background and the discontinuity of thermal evolution 
at this value are due to the presence of non-local thermodynamic equilibrium (non-LTE) to 
LTE transition of H2 ro-vibrational level population at ~ lO^cm"^. When the gas density is 
higher than this value, the cooling rate saturates and more H2 is needed to compensate for 
compressional heating. In addition, after the LTE is reached, collisional dissociation rate is 
enhanced owing to a large H2 level population in the excited levels. Thus, if a strong FUV 
radiation delays H2 formation and cooling until the critical density for LTE is reached, a 
fraction of the remaining H2 is collisionally dissociated. Thus the gas cloud is no longer able 
to cool by H2 even at a later phase of the evolution. On the other hand, if the UV background 
is slightly smaller than the threshold, the cloud begins to cool by H2 and the temperature 
begins to fall before the collisional dissociation effect becomes significant (see Figure 5b in 
Omukai 2001 for cooling rates by each process in such a case). The lower temperature allows 
further H2 formation and resultant cooling. The cooling proceeds in this accerelated fashion 
and the temperature eventually reaches the molecular cooling track. This is the origin of the 
dichotomy between the atomic and molecular cooling tracks. To summarize, the main effect 
of the FUV radiation is to photodissociate H2 directly and to decrease the H2 formation 
rate through photodissociating H~. If these two processes inhibit H2 formation and cooling 
until the critical density for LTE is reached, the gas remains warm (> several thousands K) 
and H2 is collisionally dissociated at higher densities. Thus the high density evolution is not 
affected by the presence of the FUV field and depends only on the temperature at the H2 
critical density. 
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3.2. Metallicity Effects on Irradiated Clouds 

In this section, we sliow tlie effects induced by the presence of metals and dust grains on 
the thermal evolution of gas clouds irradiated by a FUV field with a mean intensity larger 
than J2i,thr- In what follows, the total metallicity is expressed relative to the solar value, 
as [M/H] = log{Z/ZQ). Unless specified otherwise, the fractions of metals in the gas phase 
and in dust grains are assumed to be the same as in the interstellar medium (ISM) of the 
Galaxy. Specifically, the number fractions of C and O nuclei in the gas phase with respect to 
H nuclei are |/c,gas = 0.927 x 10~*Z/Zq and yo,gas = 3.568 x 10^^Z/Zq. The mass fraction 
of dust grains relative to the mass in gas is 0.939 x 10^^Z/Zq below the ice-vaporization 
temperature (Tgr < 100 K). 

In Figure [5] we present the thermal evolution of clouds with metallicity in the range 
—6 <[M/H]< —3 irradiated by extremely strong FUV radiation fields. The parameters of 
the radiation fields are (a) T, = 10^ K, J21 = 10^ and (b) T, = 10^ K, J21 = 3 x 10^ 
respectively. Under these conditions, the clouds would collapse only via atomic cooling 
in the absence of metals or dust grains (see Figs. [1] and [2]). For a metallicity as low as 
[M/H] < —6, the predicted thermal evolution follows the metal-free track. In both panels 
of Figure [5l deviations from the metal-free tracks start to appear at a density ~ lO^^cm^'^ 
when the metallicity is [M/H]~ —5.3. For the sake of comparison, thin lines show the 
expected evolution in the absence of radiation for the same initial values of metallicity. 
At metallicity [M/II]= —5.3, although the temperature drops and eventually reaches the 
molecular-cooling track at ~ lO^^cm"'^, this arrival is after the minimum in the molecular 
cooling at ~ lO^^cm"^. With a slightly higer metallicity of [M/H]= —5, this arrival takes 
place at ~ 10^^~^^cm~'^, and the temperature subsequently decreases to the minimum in the 
no-radiation case. For higher metallicities, the temperature drop occurs at lower density and 
the temperature minima becomes lower. In Figure [6] we show the cooling and heating rates 
contributed by each process during the evolution of the cloud with T* = 10^ K, J21 = 10'^ 
and [M/H] =-5. Up to 10^°cm~^, cooling is dominated by the H line emission (denoted as 
"H" in the Figure; < lO^cm"^) and H~ f-b emission ("H^ f-b"; > lO^cm^^), and the cloud 
collapses along the atomic cooling track (see Figl5]a). However, at a density ~ 10^°cm~'^, 
cooling by the dust grain ("grain") becomes dominant and causes the sudden temperature 
drop. Now the temperature is lower than that in the atomic cooling track, the H2 collisional 
dissociation rate is also reduced, which causes a high equilibrium value of the H2 fraction. 
As a result of H2 coohng, the temperature decreases further, although this effect is almost 
completely balanced by heating due to H2 formation ("H2 form"). Note that fine-structure 
line cooling ("CII, 01") is not important at such low metallicities (see the discussion below). 
Eventually, the thermal evolutionary tracks reach those of the corresponding metallicity in 
the no-radiation case (shown as thin curves in FigJS]) and evolve along them thereafter. 
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We also run models with an external UV field with parameters T,, = 10^ K, J21 = 10^, 
which is 10 times stronger than that considered in Fig|5] (a) and we found that the critical 
value of the metallicity at which deviations from the metal-free evolution appear, Zcr ~ 
5 X 10~^Zq does not depend on the intensity of the FUV radiation as long as J21 > J2i,thr- 
Furthermore, for the same value of the metallicity, the evolutionary tracks at high densities 
(> lO^cm"^) are independent of the type of external radiation, as can be seen comparing the 
results in panels (a) and (b) of Figure [51 In fact, once the evolution has reached the density 
at which molecular and atomic cooling tracks bifurcate (> lO^cm"^), coUisional processes 
rather than radiative ones dominate the energy balance (see also § I3.ip . 

It is interesting to stress that the physical processes responsible for the origin of a critical 
metallicity and its numerical value are the same as those found in the absence of FUV 
radiation (Schneider et al. 2003; Omukai et al. 2005). This is because despite the higher gas 
temperatures induced by the presence of a strong FUV field (several thousands of degrees), 
the dust temperature remains at a few tens of degrees until the energy transfer rate from gas 
to dust by collisions become important and the dust and gas temperatures approach each 
other (with the associated gas cooling). Nevertheless, the ultimate fate of protogalaxies at 
low metallicity can be significantly affected by the presence of the UV flux (see discussion in 
§ 13.41 below). In Figure [71 the evolution of dust and gas temperatures is shown for the lowest 
metallicity tracks presented in Figure[5](a). The disappearance of the dust temperature curve 
for [M/H] = —6 at a density of ~ lO^^cm"^ is due to complete vaporization of grains, which 
occurs at ~ 1300 K. At dust temperatures higher than this value, grains are no longer present 
in the cloud. Other smaller discontinuities in the dust temperature also reflect vaporization 
of some dust compounds. For example, those at ~ 130 K (at ~ lO^cm"^ for [M/H] = —6 
and -5; at ~ lO^^cm"'^ for [M/H] = —4) are due to vaporization of water ice. This figure 
indeed shows that despite the high gas temperature, the dust temperature remains low at a 
few lOK, which allows survival of grains until very high densities. 

As discussed above, 01 and CII line emission contributes negligibly to gas cooling in 
the metallicity and density range where the effects of dust grains start to become relevant 
{Zcr ~ 5 X 1O~^Z0, TT-H ~ 10^°cm~^). Metal-line cooling causes a deviation from the metal- 
free atomic track only when the metallicity reaches [M/H] ~ —3. In this case, since the 
temperature track converges to the J = track before the dust-cooling phase, two tem- 
perature minima appear at lO^cm"^ and 10^°cm~^ (see Fig|5| a, b). In the absence of dust 
grains, a higher fraction of metals is required to cool the gas at a rate such that the thermal 
evolution deviates from the atomic cooling metal-free tracks. To demonstrate this, we have 
performed a numerical experiment where we have suppressed the contribution of dust grains 
to the energy balance of the collapsing clouds. The results for models with radiation field 
parameters of {T^ = 10^ K, J21 = 10'^) are shown in Figure [H] When the metallicity is 
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below [M/H] ~ —3.5, the temperature evolution is exactly the same as the metal-free one 
(shown by the 2 x 10~^ Zq track in the Figure). For higher metallicity, fine-structure line 
cooling becomes dominant when nu < lO^cm"^ and the temperature drops abruptly by more 
than two orders of magnitude. Therefore we find that the critical metallicity [M/H] ~ —3.5 
(~ 3 X 10~^Zq) required to modify the thermal evolution is almost two orders of magnitude 
higher than in models with dust. This level of metallicity is approximately the same as the 
value at which metal cooling rate exceeds the H2 cooling rate in clouds which are already col- 
lapsing by molecular cooling (Bromm & Loeb 2003b, Santoro & ShuU 2006; Frebel, Johnson 
& Bromm 2007). 



3.3. Fragmentation Properties 

The thermal properties of star-forming clouds have an important infiuence on how they 
fragment into stars (Larson 2005). There is observational evidence that proto-stellar cores 
have a mass spectrum which resemble the stellar initial mass function (IMF), indicating that 
cloud fragmentation must be responsible for setting some fundamental properties of the star 
formation process (e.g., Motte, Andre, & Neri 1998; Lada et al. 2007; for the recent reviews, 
see Bonnell, Larson & Zinnecker 2006 and Elmegreen 2008). 

Roughly speaking, fragmentation occurs efficiently when the effective adiabatic index 
7 = dlnp/dlnp < 1, i.e., during the temperature drops, and almost stops when isothermality 
breaks (7 > 1) as also shown by the simulations of Li, Klessen & Mac Low (2003). Thus, 
consistent with Schneider et al. (2002, 2003, 2006) we can adopt the density at which the 
equation of state first becomes softer than 7 ^ 1 to identify the preferred mass scale of the 
initial fragments (Inutsuka & Miyama 1997, Jappsen et al. 2005). The fragment mass is 
given roughly by the Jeans mass (or Bonnor-Ebert mass) at this epoch, 

Mfrag = Mj(nfrag, Tfrag) oc TjJ{g n^^ll,^. (31) 

In the absence of an external FUV radiation field, the temperature of metal-free clouds 
decreases with density in the range lcm~^ ^ ""-h ^ lO^cm"^ and increases at higher densities, 
after the major coolant H2 has reached the LTE. Dense cores form around this density with 
typical masses of IO^Mq, which is close to the Bonnor-Ebert mass at this thermal state 
(Bromm, Coppi, & Larson 1999, 2002; Abel, Bryan, & Norman 2002). As the metallicity 
increases to Z^ = 10~^Zq, dust -induced fragmentation leads to solar or sub-solar fragments 
(Schneider et al. 2006), making a fundamental transition in the characteristic mass scales of 
proto-stellar cores. 



-17- 



It is important to stress that the presence of a temperature dip in the thermal evolution, 
and the softening of the equation of state 7 < 1, imply only the possibility of fragmentation. 
For example, fragmentation depends also on the initial conditions, and requires the existence 
of sufficiently large initial density perturbations. In the first cosmological objects, which are 
barely able to cool and collapse, fragmentation can be less efficient. Still, self-gravitating 
cores of mass comparable with that predicted by the above criterion are observed to form in 
high-resolution 3D simulations (Abel, Bryan, & Norman 2002, Yoshida et al. 2006). Even 
for turbulent molecular clouds of solar metallicity, 3D simulations show that fragmentation 
is efficient when 7 ^ 0.7 and it is suppressed after 7 increases to ~ 1.1 (Jappsen et al. 2005). 

The evolution with density of the effective adiabatic index, 7, is presented in Figure[9]for 
clouds with initial metallicities [M/H]=— cxd, -6, -5.3, -5, -4, and -3, irradiated by a field with 
parameters T* = lO'^K and J^ = 10^, whose temperature evolution is shown in Fig. [5] a. The 
application of the above arguments to predict the typical fragment mass from the thermal 
evolution of the clouds is not straightforward because, along the metal-free atomic cooling 
tracks and over a broad density range 10^~^^cm~^, the effective adiabatic index remains 7 ~ 
1, although slightly below unity (0.95 - 1; see Fig. [9l top panel). If we adopt 7frag = 1 as the 
threshold value of the effective adiabatic index for fragmentation, in this case fragmentation 
would be expected to occur up to densities of ~ lO^^cm"^, leading to solar-mass fragments, 
as discussed by O2001 and Omukai & Yoshii (2003). In contrast, the numerical simulations 
by BL03 show that down to the highest density reached by the simulations (< lO^cm"'^) 
fragmentation is very inefficient. Even with some degree of rotation, the cloud fragments 
at most into two pieces, resulting in a binary system. Although fragmentation might occur 
at higher densities, in BL03's calculations neither efficient fragmentation leading to the 
formation of a star cluster, nor the growth of elongation of the clouds is observed. We 
speculate that this result is due to the following reasons. The objects considered by BL03 
are those only marginally able to collapse by atomic cooling, and thus are initially close to 
the hydrostatic equilibrium. During this initial epoch, the Jeans mass is large, and density 
and velocity perturbations are erased by pressure forces. In addition to this little initial seed 
perturbation, since 7 is only slightly below unity, the growth of perturbation would be very 
slow. Thus the perturbation might not grow enough to cause fragmentation. 

Note however that, although we find that along the atomic cooling tracks H~ cooling 
is the dominant cooling agent at high densities, > lO^cm"^, this process is not considered 
in the simulation of BL03, which implements only H Lya cooling. To check whether this 
omission might cause the lack of fragmentation, we have followed the evolution of a metal- 
free cloud under the influence of an external FUV radiation fleld with T^, = lO^K and 
J21 = 10'^ but turning off the H~ cooling by hand. The result is shown in Figure [TOl With 
no H" cooling, the cloud follows a slightly higher temperature track when the density is 
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> lO^cm"^. However, below ~ lO^^cm"^, the difference is small and it has a weak effect on 
the cloud dynamics. Therefore the inclusion of H~ cooling would not affect the results of 
BLOS's simulation, which is limited to densities < lO^cm"^. 

On the basis of these considerations, we assume that for metal-free clouds irradiated by 
a strong FUV background, fragmentation does not occur during the atomic-cooling phase, 
where 7 — 1, and it occurs only when the temperature drops more rapidly, where 7 < 7fi.ag < 
1, by molecular cooling, that is when J^ < Jy^thv 

In the metal-enriched, irradiated clouds we studied, the temperature dip due to dust 
cooling occurs at very high densities, lO^'^cm"^, deep in the interior of the collapsing clouds, 
where pre-existing density perturbations might also be erased by pressure forces. However, 
in this regime we find 7 < 0.5 as shown in Figure fragmentation has also been confirmed 
to occur in two independent hydrodynamic simulations of collapsing clouds not irradiated 
by external FUV fields (Tsuribe & Omukai 2006, Clark, Glover, & Klessen 2008). Therefore, 
we expect that for metallicities [M/H] > —5, when the thermal evolutionary tracks shown 
in Figure suddenly deviate from the atomic-cooling track, in other words, when 7 falls 
sufficiently below unity (Figure [H]), the clouds begin a vigorous fragmentation, which then 
lasts until the temperature increases again. The value of 7frag to cause fragmentation is 
uncertain as discussed above, but likely to be slightly below unity. In the following, for the 
sake of definiteness, adopt 7frag = 0.8 as the fiducial value below which fragmentation is 
triggered (the lower horizontal lines in Fig. [H]), and use it to define the properties of the 
fragments. This choice refrects the fact that for 7 = 0.7 efficient fragmentation has been 
observed to occur in the numerical simulations of Jappsen et al. (2005). In all cases, once 
molecular cooling becomes efficient, 7 soon falls below < 0.5 (Fig. [9]). Varying the threshold 
value 7frag, say, by ~ 0.1, leads to fragmentation densities whose differences are within an 
order of magnitude. For example, when [M/H]= —4 fragmentation begins at log nH(cm~^) = 
7.5, 8, and 8.3 for 7frag = 0.7, 0.8, and 0.9, respectively. We assume that fragmentation stops 
when 7 exceeds unity again. For [M/H]= —4, this occurs at tih ~ lO^^cm"^. 

The mass scale of the final fragments is given by the Jeans mass at the temperature 
minimum, i.e., when 7 exceeds unity. When the initial metallicity is [M/H] ^ —5, the tem- 
perature minimum corresponds to 300 K at wh = lO^^cm"^, and thus the typical fragment 
mass is O.IM©. As the metallicity increases, both the density and temperature at the frag- 
mentation scale decrease, being (lO^^cm"^, 150K) for [M/H] ~ —4, and (lO^^cm"'^, 30K) for 
[M/H] ~ —3. However, the corresponding fragment mass scale remains ~ O.IMq, because 
the variations of density and temperature almost cancel out (see eq. EI]). In some cases, e.g. 
[M/H] ~ —3 gas in a J21 > t/2i,thr field, two fragmentation epochs (log nn = 3.5 — 5.1 and 
8.5 — 10.7 for [M/H] = —3) appear, which corresponds to two dips in the temperature (or 
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7) evolutionary track. The outcome of this kind of track is not clear without any numerical 
work studying their effect. Here, we speculate that the first dip produces clumps as a result 
of the fragmentation of clouds. Then the clumps fragments again into cores owing to the 
second dip. 

In the absence of dust, the temperature minimum appears at a lower density, nn ~ 
lO^cm"^, and higher metallicity [M/H] ^ —3.5 (see Fig. [8]). Therefore, the corresponding 
fragment mass remains as high as 10 — IOOMq and the formation of sub-solar mass fragments 
is not possible in this case. This property of pre-stellar clouds enriched only by gas-phase 
metals has been already proven to hold in the absence of external FUV fields (Schneider et 
al. 2006). 

To summarize, our results show that in the presence of a sufficiently strong FUV radia- 
tion field the collapse of metal-free clouds by molecular cooling is inhibited and it can proceed 
only via atomic cooling. Under these conditions, cloud fragmentation is highly inefficient, 
leading at most to the formation of a binary system. The typical mass of pre-stellar clouds 
is therefore 10^~^Mq and the formation of a super massive star, seed of a super massive 
black hole, is the likely outcome of the evolution (BL03). However, this scenario is altered 
as soon as trace amounts of metals and dust grains are present in the collapsing clouds: dust 
cooling leads to fragmentation of the clouds into sub-clumps with mass as low as ~ O.IM© 
already at a floor metallicity of Zcr ~ 5 x 10~^Zq. This conclusion holds independently of 
the intensity and spectrum of the FUV radiation fleld. In the absence of dust, an enrichment 
level of Zcr ~ 3 X 10~^Zq is required for 01 and CII line cooling to fragment the cloud; the 
fragments in this case are predicted to be relatively more massive, ~ 10 — lOOM©. 



3.4. Dynamical Interactions and Accretion 

Since dust-induced fragmentation takes place at high densities, a dense proto-stellar 
cluster is expected to form (Omukai et al. 2005, Schneider et al. 2006, Clark et al. 2008). As 
an example, when the initial metallicity of the collapsing cloud is [M/H] = —5, the sudden 
temperature drop, where 7 < 7frag = 0.8, begins at Tdrop ~ 3500 K and ndrop — lO^^'^cm"^. 
At this stage, the size and mass of the cooling region, or proto-cluster, are given by the 
corresponding Jeans length, Aj ~ 4 x 10^^ pc, and mass M^i — 70Mq. When a different 



threshold value 7fi.ag is adopted, these quantities change, e.g., to Aj ^ 3 x 10 pc, and M^ 



cl 



4OM0 for 7fi.ag = 0.7 (0.9, respectively). In the following order of magnitude estimation, we 
use Aj ~ 1 X 10^^ pc and Md ~ IOOMq as typical values. After virialization, the proto- 
stellar cluster has a size half of this. Since each ultimate fragment has a typical mass of 
^frag ~ O.IM0, which is set by the Jeans mass at the end of the fragmentation process 
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(nn ~ 10"^^cm~^), we expect that up to A^* ~ Md/Mfrag ~ 1000 low-mass star can be formed 
and confined into a small region of size ~ 0.01 pc. The difference between the formation 
epochs of each protostar is of the order of the free-fall time of the proto-cluster gas. Since 
the cluster begins to form in a dense cloud with density ~ 10^°cm~^, protostar formation is 
synchronized on a timescale of ~ 300 yrs. 

The fate of dense, compact star clusters has been discussed extensively in the literature 
(see, e.g. Rasio et al. 2004 for a recent review, focusing on the possibility of intermediate BH, 
IMBH, formation through a runaway collapse that is relevant in our case). It is important 
to stress that, even assuming a star formation efficiency of order unity (which seems likely 
when the density exceeds > 10*^cm~^; Alves, Lombardi & Lada 2007), the stellar IMF will 
be strongly affected by gravitational interactions, collisions and mergers. In fact, observed 
properties of present-day star forming regions, as well as numerical simulations, suggest 
that gravitational fragmentation is probably responsible for setting a characteristic stellar 
mass but the full mass-spectrum and the Salpeter-like slope of the IMF are most likely 
formed through continued accretion and dynamical interactions in a clustered environment 
(see Bonnell et al. 2006 and references therein). Furthermore, in young and compact star 
clusters supermassive stars may form through repeated collisions (e.g. Portegies Zwart et al. 
1999, 2004; Ebisuzaki et al. 2001). 

We can therefore ask, what is the expected fate of the dense star cluster forming in our 
clouds? The evolution of a star cluster with half-mass radius i?ci and mass M^y proceeds on 
the dynamical friction timescale (Binney & Tremaine 1987), 

1.2 X 10^ / r Vf i?ci Y^'"^ f Mci V^ f m. 



^'"^- InA ^'Vo.OlpcJ VO.OlpcJ yi^mj \lMrJ ' ^^^^ 

which is the time required for a star with mass m*, which is on the massive side of the 
spectrum, to sink from the radius r to the cluster center by gravitational interactions with 
background, less massive stars. In the following, the Coulomb logarithm InA is taken to be 
~ 7, a value typical for open clusters. In the above equation, the density profile is assumed 
to be isothermal, p oc 1/r^. The dynamical friction timescale was originally derived for a 
fixed background. What actually occurs for a cluster on this timescale is the equipartition 
of kinetic energy among the member stars. Heavy stars move slowly and then drop deeper 
in the potential well, leading to mass segregation in the cluster. The stellar merger rate is 
greatly enhanced if higher-mass stars reach the cluster center within the lifetime of a very 
massive star, i.e. if tfnc is less than a few Myr. 

Using the estimated size and mass of the proto-cluster at the onset of dust-induced 
fragmentation for the [M/H] = -5 track, the dynamical friction timescale for a star initially 
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at radius r can be rewritten more generically as 



where M^ is the mass enclosed inside radius r and we have expressed Md and Rc\ in terms 
of ndrop and Tdrop (note that ridrop drops out of the equation). Following Portegies Zwart et 
al. (2004), we assume that a very massive star can be formed by stellar mergers if, 

tMe<4Myr. (34) 

From equation (l33l) , we infer that the inner region of mass 

Mmc = 5 X lO^Me {^) '"' f ^^1 ''' (^ ''' (35) 

°V1.22/ V5000Ky \IMq) ^ ' 



satisfies the condition fl34j) . The mass fraction of sinking stars relative to the cluster back- 
ground stars, /sink, is uncertain and probably depends on the stellar mass spectrum, but 
can be safely assumed to be less than a half. In addition, not all the stars that sink to 
the center are incorporated into the runaway merging object. Since merging events usually 
proceed via three- or more body interactions, where the runaway object and a star coalesce 
by kicking the lightest star (Portegies Zwart et al. 1999), the merging efficiency, /merg, with 
the runaway object among stars fallen to the cluster center can be assumed to be about a 
half. Thus, the mass of the central object resulting from this process can be estimated as, 

-'"cen ^ /sink/merg-'"fric \'~'^) 

^ 3.5 X IO^Mq (^-^ (^^^^ fTAro^\ f^!^\ .3-^) 

° Vo-25y V 0.5 y V5oooKy \lMrJ ^ ' 

As can be seen in Figures [5] and [9l dust-induced fragmentation in clouds with a metal- 
licity [M/H] < —4, starts at Tdrop ~ 5000 K, ridrop ^ lO'^cm"^ and the corresponding 
proto-cluster mass is Md < 1O^M0. That is, Md < Mmc and the entire cluster satisfies the 
condition fl5^ . In this case, the mass of the central object is limited by the mass of the cluster 
rather than by Mfnc and its final mass can be estimated as, Mcen = /sink /merged ^ IOOMq. 
On the other hand, dust-induced fragmentation of clouds with metallicity, —4 < [M/H] < 
—3, leads to proto-stellar clusters with masses M^x > M^^-^c and the mass of the central 
object is given by equation (I57|) . Thus, in this case a very massive star can form with 
M:en ^ 35OM0. In either case, we note that the metal-poor, massive star ultimately form- 
ing at the center of the halo is likely to leave behind a seed BH remnant - except in a narrow 
range of metallicity, where they produce a pair instability supernova - either by direct col- 
lapse or by fallback (Heger et al. 2003). For higher metallicity, there are two episodes of 



-22- 



fragmentation; a metal-induced one at low-density and a dust-induced one at high-density. 
In this case, the mass scale of the star cluster is relatively low and a massive BH seed is not 
formed. 

As pointed out above, the critical metallicity levels we find in the case of strong FUV 
irradiation, at which the cloud behavior is modified from the metal-free case, is very similar 
to the critical metallicity found in earlier work for J = 0. It is therefore important to ask 
whether the presence of the flux will, in fact, make any difference to the ultimate fate of 
the cloud. At metallicities above [M/H] ~ —3, the flux has essentially no impact on the 
evolutionary track of clouds at high densities. For example, comparing the thin and thick 
solid curves of Figure [5] (a) it is clear that when [M/H] = —3, the flux has no effect at 
nn > lO^cm"^. At lower densities, we expect that the first fragmentation phase will occur 
as the Jeans mass drops from Mj ~ 10^ M0 to ~ 10^ M© at nn ~ lO^cm"'^ in the J = 
case, and to ~ 10^ M© at lO^cm"'^ in the strong UV background case. Therefore, the size of 
the molecular clumps that form is larger in the J = case. However, the thermal evolution 
thereafter is the same: the molecular clumps will experience a second phase of fragmentation 
when the Jeans mass falls further, from Mj ~ 10 M© to ~ 0.1 Mq. Thus, for protostellar gas 
clouds with [M/H] ~ —3, the presence of an external UV field determines only the amount 
of gas in the envelope in which the ~ 10 Mq star cluster is embedded. 

In contrast, the presence of a FUV background significantly affects the evolution of 
protostellar clouds with lower values of metallicity, [M/H] < — 3. In these models, the 
atomic gas experiences only the second phase of fragmentation induced by dust cooling, and 
therefore more massive star clusters (10^ — lO^M©, increasing with metallicity) are formed 
compared to those in the J = limit (a few - IOMq, depending weakly on metallicity). In 
addition to the size of the star clusters, the mass and physical conditions of the corresponding 
envelopes are different: when J = 0, the star cluster is embedded in a molecular envelope 
of 10^~^ Mq with temperature of a few 100 K, while, in the presence of an external UV 
field, the surrounding envelope is more massive (10^~® M©) and it is made by atomic gas at 
several 1000 K. As a consequence, the presence of a strong UV background case favors the 
formation of a more massive central star by stellar merger (larger stellar cluster mass and 
higher Tdrop in eq. EZ])- 

Finally, our results and the discussion above suggest that the direct formation of a 
supermassive star or SMBH as massive as ~ 10^"^ M0, as envisioned in the metal-free case, 
is not possible when the metallicity is above a critical value, and the gas fragments into 
smaller pieces. However, we note that if fragmentation of the inner regions of the collapsing 
protogalaxy is not fully efficient, and if radiative and mechanical feedback from the stars 
does not expel the leftover gas from the nucleus, then the star cluster at the center of 
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the halo, and its coalesced massive remnant star, can be embedded within a thick residual 
gaseous envelope with temperature Tgnv ~ several thousand K. Since this gas envelope is 
self-gravitating, with a temperature below the virial temperature, it can undergo dynamical 
collapse, and may still produce a SMBH either directly or by accretion onto the central 
stellar-mass BH at the Bondi rate (see, e.g. Begelman et al. 2006 and discussion therein). 
In particular, in the UV irradiated case, the accretion rate onto the star cluster, and possibly 
to the central coalesced star, is high - this is because of the high temperature in the envelope, 
and since the accretion rate of the self-gravitating gas is given by (Shu 1977): 

M ~ -^ (38) 

Gr 

(rp \ 3/2 

-^^^ . (39) 

soook; ^ ' 

In conclusion, higher initial mass, higher accretion rate, and larger amount of reservoir gas 
are more favorable for BH growth in the strong UV case than the case without radiation. 



4. Summary and Conclusions 

In this paper, we have investigated the thermal evolution and fate of proto-stellar gas 
clouds in Tvir > 10^ K halos irradiated by a strong FUV background. Under these conditions, 
which may apply to some dwarf galaxies collapsing close to the epoch of reionization, we 
find that: 

• The effect of an external UV background is to photodissociate H2 directly and to 
decrease the H2 formation rate through photodissociation of H~. When the UV back- 
ground reaches a critical threshold value, J2i,thr, these two processes inhibit II2 forma- 
tion and cooling until the critical density for LTE is reached. Thereafter, the gas cloud 
can cool only via atomic hydrogen transitions. 

• For gas clouds of primordial composition, an external UV background with intensity 
J21 < J2i,thr only delays the onset of II2 formation and H2 cooling becomes important 
at some (higher) density: fragmentation occurs at densities lO^cm"^ — lO^cm^^ leading 
to average fragment masses in the range IO^Mq — lO^M©, similarly to the case with 
J21 = 0. 

• For J21 > J2i,thr, not enough H2 is formed to activate cooling and the evolution of 
primordial clouds is controlled by atomic (H and H^) cooling. The clouds collapse 
nearly isothermally with a temperature of several thousands K ("atomic track") up 
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to very high densities ~ lO^^cm^^. According to previous numerical calculation by 
BL03, the clouds collapse directly into a single 10^ — IO^Mq object, leading to super 
massive star and SMBH formation. A core-envelope structure inevitably develops 
under these circumstances, and the star grows by accretion onto an initially small 
inner core. During the accretion phase, radiative and mechanical feedback effects might 
become important and eventually halt the accretion at some phase (e.g., Omukai & 
Palla 2003; McKee & Tan 2007). If so, the mass of the central object can remain far 
©• 



below lO^-^MfT). 



• Independently of the values and properties of the external FUV field (as long as it is 
J21 > J2i,thr)5 deviations from the metal-free "atomic track" start to appear when the 
gas is enriched by even trace amounts of metals and dust. When Z > Z^r ~ 5 x 10~^ Zq, 
dust cooling induces fragmentation at nn ~ lO^'^cm"'^ and a proto-stellar cluster is 
expected to form with average proto-stellar (fragment) mass of ~ O.IM0. If only 
gas-phase metals are present, a two orders of magnitude larger value of metallicity is 
needed, Z^r ~ 3 x IO'^^Zq, before CII and 01 line-cooling induce a deviation from 
the "atomic track", leading to fragmentation at nn ~ lO^cm"^ and to proto-stellar 
clusters with average proto-stellar (fragment) mass of 10 — IOOMq. 

• The physical processes responsible for the origin of a critical metallicity and of its 
numerical value are the same as those found in the absence of an external FUV field. 
However, due to the higher gas temperature, the final outcome of the proto-stellar 
cloud collapse can be significantly affected. Namely, if we assume that the size of the 
proto-stellar cluster formed by dust-induced fragmentation is set by the Jeans mass 
at the onset of the rapid temperature drop, it depends on the intensity of the FUV 
background field: when J21 < J2i,thr, relatively small clusters are formed (a few - IOM0) 
whereas when J21 > J2i,thr, very dense star clusters with masses 100 — 1000 Mq are 
formed, at the center of which stellar coalescences are expected to occur. The central 
merger object might grow to a very massive star of a few 100 Mq. 

• In addition, the presence of an external FUV background affects the physical conditions 
of the envelope surrounding the proto-stellar clusters: when J21 < J2i,thr the envelope 
is fully molecular, with a mass of 10^ — 10'^ Mq and a temperature of a few 100 K. 
Conversely, when J21 > J2i,thr the envelope is made of atomic gas and reaches a mass 
of 10^ — 10^ Mq and temperature of several 1000 K. The higher temperature and larger 
gas reservoir favors BH growth by accretion, which can be as high as lO^^MQyr^^. 

According to the above, the conditions that would allow the formation of the direct 
formation of a SMBH are (i) to be hosted within a T^w ~ 10"^ K halo, which is ii) irradiated 
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by a strong UV field with J21 > t/2i.thr, and (iii) still metal and dust free, with Z < Zcr- The 
main new result of the present paper is that if the metalicity is too high, so that condition 
(iii) does not hold, then instead of a SMBH, a dense cluster of low-mass star forms at the 
halo nucleus. The stars in such a cluster may still rapidly coalesce into a single massive 
star, which may produce an intermediate-mass BH remnant, but with a smaller mass of 
M< 10^- 10^ Mq. 

While the above conclusion that even trace amounts of dust enable cooling and frag- 
mentation of the proto-stellar clouds appear to be robust, the exact value of the threshold 
metallicity is vulnerable to the uncertain nature of dust in early protogalaxies (see also 
Schneider et al. 2006). However, it is interesting to note that when Z > Z^r the characteris- 
tic fragment mass - which is related to the characteristic stellar mass - is highly insensitive 
to environmental conditions, such as the presence of an external FUV radiation field, as also 
recently discussed by Elmegreen, Klessen & Wilson (2008). 

We warn the reader that our discussion on the nature and evolution of the resulting 
proto-stellar cluster is still speculative as it is based on a few numerical experiments which 
apply to dense stellar systems in present-day star forming regions (see the discussion and 
references in section [37^ . For example, we assume that the size of the cluster is set at the 
onset of the efficient cooling phase, which appears plausible but has not yet been confirmed. 
In particular, in models with J21 > J2i,thr and Z > 10~^Zq (when both dust grains and 
gas-phase metals are present), the thermal evolution curves appear to have two separate 
temperature minima, which correspond to metal- and dust- induced cooling and fragmen- 
tation. The fate of these collapsing clouds is at present unknown, and dedicated numerical 
simulations would be highly desirable. If proto-stellar clusters are indeed formed under the 
conditions that we suggest, the formation and the nature of a central object by repeated 
collisions and accretion is highly uncertain. For example, the fraction of stars on the mas- 
sive side of the spectrum which falls to the cluster center (i.e., /sink in eq. [2ZD by dynamical 
friction is unknown and may vary significantly depending on the IMF. 

Despite the above uncertainties, our results suggest that even trace amount of metals 
preclude the rapid formation of SMBHs as massive as M ^ 10^^^ Mq in protogalactic halos. 
While such promptly appearing SMBHs would help solving the puzzle of the M > 10^ Mq 
quasar black holes at z > 6, our results suggest that the low-metallicity halos may instead 
produce dense stellar clusters - the cluster may coalesce to produce an IMBH, but still with 
a much lower mass oi M ^ 10^~^ Mq. 

This study is supported in part by the Grants-in-Aid by the Ministry of Education, 
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number density log n^^ (cm"^) 



Fig. 1. — Temperature evolution of metal-free clouds irradiated by a UV flux. The spec- 
tral shape is that of a black-body spectrum with T* = lO^K. Models are shown with 
FUV intensities at the Lyman limit of J21 = 0, 1, 10, 100 and 10^, in the usual units of 
10~^^erg cm~^ sr~^ s~^ Hz~^ (solid and dashed curvesfrom bottom to top; see the legend 
in the panel). Diagonal dotted lines correspond to different constant Jeans mass. Models 
with higher initial temperature (210 K in contrast to 21 K in the fiducial models) are also 
shown by dotted lines. For those with J21 = 10, 100 and 10^, the temperature evolution at 
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nn ^ 10^'^cm ^ completely overlaps with that predicted by the fiducial models. 




number density log n^^ (cm ^) 



Fig. 2. — The same as Figure [H except assuming a harder spectrum, with T^ = lO^K and 
intensities J21 = 0, 10, 100, 10^, 10^, 10^, and 3 x 10^ (sohd and dashed curves from bottom 
to top). 
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Fig. 3. — Contributions of various processes to the total cooling rate, as a function of 
the number density, for a metal-free cloud irradiated by an extremely intense FUV field 
(T* = lO^K, J21 = 10^; with the temperature evolution shown in Figured]). The meaning of 
the symbols is as follows: "compr" indicates compressional heating; "H~ f-b" cooling by H~ 
free-bound emission; "H" cooling by H line emission; "H2" cooling by H2 line emission. 
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Fig. 4. — H2 and H~ photodissociation rate coefficients as a function of radiation temperature 
for J21 = 1. Also sliown for reference are tlie Habing parameter Gq and dilution factor W. For 
a fixed intensity J21 = 1 at 13.6 eV, the H2 photodissociation coefficient is almost constant 
with T* while the H~ coefficient is about four orders of magnitude higher for T^, = lO^K, and 
thus reduces H2 formation by a large factor (see text). 
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Fig. 5. — The temperature evolution of clouds with initial metallicity [M/H]= -6 (solid), 
-5.3 (dotted), -5 (short-dashed), -4 (long-dashed), and -3 (dash-dotted) irradiated by a FUV 
field with (a) T, = lO^K, and J21 = 10^ and (b) T, = lO^K, and J21 = 3 x 10^ Thin curves 
show the results obtained without an external FUV field for initial metallicities [M/H]= -5.3, 
-5, -4, and -3 (the same line types as the irradiated cases with the same metallicity). Due 
to photoelectric heating, when [M/H]=-3, the temperature at the lowest densities is higher 
than in the other models. 
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Fig. 6. — Cooling and heating rates contributed by each process during the collapse of a cloud 
with a metallicity of [M/H]=-5 and external FUV radiation with (T* = lO^K, J21 = 10'^). 
The corresponding temperature drop at rin ~ 10^°cni~^ (shown in Figure E]) is caused by 
dust cooling. The meaning of the symbols is as follows: "compr" indicates compressional 
heating; "grain" cooling by dust thermal emission; "H~ f-b" cooling by the H~ free-bound 
emission; "CII, 01" cooling by the CII and 01 fine-structure line emission; "H" cooling by 
the H line emission; "H2" cooling by the H2 line emission; "H2 form" heating by the H2 
formation. 
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Fig. 7. — The dust temperature of clouds (solid curves) with metallicity [M/H]=-6, -5, 
and -4 and with an external radiation of (X,, = lO^K, J21 = 10^) as a function of density. 
For the same models the dotted lines indicate the corresponding gas temperatures. The 
dust temperature curve disappears at about wh ~ lO^^cm"^ for [M/H]=-6, reflecting the 
vaporization of the grains. 
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Fig. 8. — Temperature evolution of clouds when all metals are assumed to be in the gas 
phase (i.e. no dust). The radiation parameters are T^, = lO^K, J21 = 10^. Models with 
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initial metallicities Z = 2 x 10 ^Z© (solid), 3 x 10 ^Z© (dashed), and 10 ^Zq (solid) are 
shown (from top to bottom). 
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Fig. 9. — Evolution of effective adiabatic indices 7 = dlnp/dlnp for the models shown in 
Figure [S](a), i.e., clouds irradiated by a field with parameters T^, = lO^K and J* = 10^, 
with metallicities [M/H]= —00, -6, -5.3, -5, -4, and -3 (from top to bottom). The cases 
of [M/H]= —00 and -6 are identical and shown in the same panel (top). The horizontal 
lines indicate 7 = 0.8 and 1. The former is adopted as the fiducial value below which 
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fragmentation is triggered. We assume that fragmentation stops when 7 becomes > 1. 
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Fig. 10. — Effects of the H~ free-bound coohng on the temperature evolution. The evolution 
of a metal-free gas irradiated by a FUV field with (T* = lO'^K, J21 = 10'^) is shown with and 
without the H~ free-bound cooling. 



